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O^l ' Chandler wobble excitation and damping, one of the open problems in geophysics, is treated 

, as a consequence in part of the interaction between Earth and a hypothetical oblate ellipsoid made 

I of dark matter. The physical and geometrical parameters of such an ellipsoid and the interacting 

torque strength is calculated in such a way to reproduce the Chandler wobble component of the 
polar motion in several epochs, available in the literature. It is also examined the consequences 
upon the geomagnetic field dynamo and generation of heat in the Earth outer core. 
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li^; I Introduction 

I A 305-day Earth free precession was predicted by Euler in the 18th century and was sought 

\^ ' since then by astronomers in the form of small latitude variations. F. Kiistner in 1888 and 

. S. C. Chandler in 1891 detected such a motion but with a period of approximately 435 

*^ I days. Later it was shown by Love [Ij and Larmor [2j that such a disagreement between 

• ■ theory and observation is in part consequence of elastic deformation of the Earth which 

. was not considered by the Euler rigid body model (see Smith and Dahlen [3j for further 

OO ! corrections explanation). The so called Chandler wobble (CW) is actually one of the major 

' components of the Earth polar motion, together with the annual wobble (AW) which has 

period of nearly 365 days. The composition of the two wobbles results in a motion with 
^ , period of around 6.2 years due to the beat phenomenon. The ever increasing measurement 

Vh ' precision, for which techniques like very long baseline interferometry (VLBI), global posi- 

tioning system (GPS) and lunar laser ranging (LLR) are employed today, has lead to the 
conclusion that both wobbles have varying amplitudes. For AW, this may be attributed 
to seasonal displacement of atmospheric and water masses. On the other hand, excita- 
tion and damping of the Chandler wobble has become a puzzle to investigators. Several 
mechanisms have been proposed for its excitation like snow, hidrological, atmospheric and 
oceanic mass displacements, and large earthquakes, with defenders and opponents arguing 
for and against repeatedly and, in view of the large energy involved, the problem is still an 
open question. In more recent papers. Gross [4j attributes the excitation to ocean bottom 
pressure variations and changes in ocean currents due to winds while Seitz et al. [5] focus 
on combined effects due to atmospheric and oceanic causes. 

Since it is believed that dark matter (DM) and dark energy composes around 96% of 
the mass of the Universe we investigate in this paper the possibility of DM to contribute, 
at least in part, for the excitation and damping of Chandler wobble. The effect of DM 
upon regular matter has been considered mainly on galactic-scale although more recently 
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some authors have studied smaUer scale effects. For instance, Proggatt and Nielsen [6]) 
have investigated the possibility of star-scale effects due to small DM balls; Frere et al. [7] 
studied the presence of a DM halo centered in the Sun and its influence upon the motion of 
planets; numerical simulations due to Diemand et al. [8j show the presence of DM clumps 
surviving near the solar circle; Adler [9] suggests that internal heat production in Nep- 
tune and hot-Jupiter exoplanets is due to accretion of planet-bound, not self-annihilating 
DM. Although the nature of the particles that form DM is not known precisely (see for 
instance Khalil and Muhoz [10]) we assume that, as usual, those particles do not interact 
with regular matter through electromagnetic force but only gravitationally. Furthermore, 
we speculate that they can interact with each other in such a way to set up an oblate 
(due to rotation) ellipsoid that might be present inside Earth, occupying the same space 
without violating any physical law. It can be shown that such an ellipsoid, made of self- 
interacting DM (otherwise such a structure would not be possible), would be attached 
to Earth through a spring-like restoring force. In this line, excitation and damping of 
CW could be, at least partially, consequence of energy exchange between the two bodies, 
associated to the restoring torque between them. The model does not consider the other 
known sources for the process nor the Earth internal structure, what are of course over 
simplifications. However this becomes possible a simple starting point for the calculations, 
to be sophisticated in a later step. Hopfner pj^ tl2j has managed to isolate the CW and 
AW components from available polar motion data. Our goal is to reproduce Hopfner's 
results for Chandler wobble in various epochs by adjusting the unknown physical and geo- 
metrical parameters of the dark matter ellipsoid. For this purpose, we describe both body 
motion by solving numerically the set of differential equations that emerge from Euler 
equations and demonstrate that CW can be reproduced in such manner. This is explained 
in the next section while results and discussion are presented in section III. 



II The model 

We start by calculating the torque that arises when the polar principal axes 3 and 3* of 
two oblate ellipsoids with isotropic densities interacting through the gravitational force 
are tilted by an angle a (see Fig. 1). The gravitational potential at a point outside the 
internal ellipsoid (see e. g. Stacey [13]) is given by 

CM* C 

V[r,e)^-^^ + ^{C* -A*)P2{cose) (1) 

neglecting higher order terms. Here M*, A* and C* are the mass, equatorial and polar 
momenta of inertia of the internal ellipsoid, respectively, and P2{cos6) = (3cos^ 6 — l)/2 
is the Legendre polynomial of second order. The torque over an infinitesimal mass dm 
located at that point is then 

dT = -dm^ = 3g(g* - A*) ^-^^26) dm. (2) 

and is directed orthogonal to the plane formed by the polar principal axis 3* and the 
position r of dm (see Fig. 1). Considering the Earth tilted by an angle a around the 
1-axis (Fig. 1) with respect to the internal ellipsoid it suffers therefore a torque given by 
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the integral 

sin (j) dr (3) 



directed towards the identical axes 1 and 1*, as projected by the presence of sm<j) in 
this integral. The integration should in principle be performed over the entire Earth 
volume. However considering the Earth density as isotropic, as it is done in the Preliminary 
Reference Earth Model - PREM (Dziewonski and Anderson [Hj), it can be shown that 
the spherical volume internal to the Earth with radius the same as the polar radius does 
not contribute to the integral. Therefore the integration over r ranges from Rp, the polar 
radius, to the radius of the geoid at a given direction. On the other hand, it is well known 
(see, for instance, Stacey [13]) that the equation of the geoid is written in the first order 
approximation in terms of the co-latitude 9' as 

rg^ail-fcos^e'), (4) 

where a is the equatorial radius and / = 1 — Rp/a is the flattening. Let us write above 
equation for the case where the geoid is tilted by an angle a around the 1-axis (Fig. 1), 
which yields 

Vg ^ a 1 — /(cosacos^ — sinasin^sinc/))^ , (5) 

where cj), together with co-latitude 6, is part of the spherical coordinates in the 1*2*3* 
reference frame. Therefore the torque given by equation (3) becomes 

r = -pG{C* - A*) r de sm{29) sine 
2 JO 

X / dcf) sm(f)ln < ^ 1 — /(cos a cos 6* — sin a sin sin (/))^ >, (6) 
JO \ 



where p is the average Earth mass density in that region, taken as 2600 kg/m^ 

Considering that / = 1/298.257 << 1, we can use the approximation ln(l + x) ~ x, 
obtaining 

r = ci sin(2a), (7) 

with ci = ^pG{C* — A*)f. Notice that this is a restoring torque when a is around 90°, 
which gives the stable angular position if no initial angular momenta were involved. 
We describe the Earth-ellipsoid motion by the Euler equations 

Adji + (C - A)uj2UJ3 = n (8) 

Adj2 - (C - A)UJIUJ3 = T2 (9) 

CuJs = T3 = =^ UJ3 = constant (10) 

A*OjI* + {C* - A*)UJ*2*UJ*3* = T* (11) 

A*dj^, - {C* - A*)ujI,uj*s, = T2* (12) 
C*uj^* = t| = =^ (^3* = constant (13) 

where A and C are the equatorial and polar momentum of inertia of Earth, which is 
supposed to be axially symmetric (B = A), uji, ijj2, uJs are the components of the Earth 
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angular velocity u in the 123 reference frame, which is attached to Earth such that the 3- 
axis coincides with the polar principal axis, i. e., the symmetry axis of the oblate ellipsoid, 
1- and 2-axis are equatorial principal axes and ri , T2 are components of the torque suffered 
by Earth. Similar internal ellipsoid quantities are denoted by an asterisk. Notice that ra 
and r| vanish, what implies that and are constant. Notice also that we are not 
considering external torques, mainly exerted by the Sun and the Moon, which produce 
the precession of the equinoxes and is of no interest here. On the other hand, when 
Ti = r2 = r * = r| = equations (8)-(13) lead to the well known Eulerian free precession 
motion. 

We write the angular velocities in terms of the Euler angles 9, (p, ip and their derivatives 
in time, defined with respect to a fixed reference frame with axes X, Y, and Z (see Fig. 2). 
This frame is constructed such that the total angular momentum of the Earth-ellipsoid 
system, which is a constant of the motion, coincides with the Z-axis. The relations are 
the following 

LOl = 9 COS ip + (f) sin 9 sin ^ (14) 
UJ2 = — ^?sinV' + (/'sin^JcosV' (15) 
uJs = (f) cos 9 + ^ = constant ^ ip = uj^ — cos 9. (16) 

Similar relations for cj|, , 0^2* , 1^3* can be written in terms of the Euler angles 9*, (f)*, ip* , 
referred also with respect to the XYZ fixed system, and their derivatives in time. Notice 
that equation (16) serves for the purpose of eliminating tp, while the same elimination can 
be done for ip*. The components of the torque r over the Earth can also be expressed in 
terms of the Euler angles 

Ti = 2ci cosa [— sin^cosV'COS^* + (cos0cos(/)cosV' — sin(;/!)sin'0)sin0*cos(/>* 

+ (cos 9 sin (p cos ^p + cos (p sin tp) sin 9* sin (p*] (17) 

T2 = 2ci COS a[sin^ sin-i/; cos — (cos ^? COS (/)sin^ + sin i}!)COS '0) sin6* COS ^* 

+ (— cos sin ^ sin '0 + cos cos tp) sin 9* sin (p*]. (18) 

The expressions for the components r* and r| of the torque t* over the internal ellipsoid 
are like above, provided the exchanges 9 -f-^ 9* , (p ^ (p* , tp tp* are made. In the same 
way, we also need the expression for cos a 

cos a = sin ^ sin 9* cos{(p — (p*) + cos 9 cos 9*. (19) 

After substituting equations (14)-(18) into equations (8), (9), (11), (12) and considering 
that we are dealing with real quantities we can eliminate tp, tp*, tp, tp*, getting the following 
set of four second-order differential equations 

C ci 
sin 6* - -i039 + 29<j) cos 9 - 2-f cos a sin 9* sin((p* - (p) = (20) 

9 - 4? sin 9cos9 + ^uj^^ sin 9 + 2^ cos a [sin 9 cos 9* - cos 9 sin 9* cos{(p - cp*)] = (21) 

^*sinr- ^w^J* + 2r0*cosr-2^cosasin0sin((/)-</)*) =0 (22) 

9* - (j)*^ sin 9* cos 9* + ^a;t ^* sin 0* + 2-^ cos a [sin 9* cos 9 

- cos 9* sin 9 cos{(p - cp*)] = 0. (23) 
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As a test for these equations we consider the case where ci = 0, what means that the 
ellipsoids are not interacting and therefore free precession occurs. It is known (see e. g. 
Symon [15]) in this case that for each ellipsoid both polar principal symmetry axis (3- 
axis) and the rotation axis, which contains the angular velocity vector cj, have precession 
around the conserved angular momentum vector L forming the so called body cone and 
space cone, with semi-angles at and a^, respectively (see Fig. 3). Notice that the 3-axis, 
u and L are always in the same plane. The constant angular velocity of such a precession 
around fixed L is 

(24) 



smob 



sma. 



A cos 6- 



7^3 



with 



(25) 



being the constant angle that gives the orientation of the principal 3-axis with respect to 
fixed L, and 



tana;, 



cos as 



W3 

1 + j3 cos^ Ofo 
l + (2 + /3)/?cos2 Ob 



/3 



~A 



1. 



(26) 
(27) 

(28) 



However, with two ellipsoids, we have a more general situation where the fixed reference 
frame Z-axis is taken to be in the direction of the total angular momentum Lj, which is 
conserved if external torques are negligible, and not of the Earth angular momentum L. 
Therefore we define a reference frame xyz, such that z is in the direction of L, oriented with 
respect to XYZ through the Euler angles Sg, 6^, 6^ (see Fig. 4). While the principal 3-axis 
has precession around z-axis with constant angular velocity and angular amplitude 6', 
the 3-axis has motion described in terms of the Euler angles 9, (j) with respect to the frame 
XYZ as 



cos ( 



tan ( 



sm ( 



sm{(pt + e) sin 5g sin 6^ — cos((/)t + e) sin do cos 5, 



+ cos 6' cos df) 



sin 6' 


Bi sm(<pt + 


e) + B2 cos{<j)t + e) 


+ cos 9' sin 6g sin 6^ 


sin 9' 


Bs sm{4>t + 


e) + B4Cos{(j)t + e) 


— cos 6' sin 5g cos 5^ 



(29) 
(30) 



with e being an arbitrary phase and 



Bi 
B2 
B3 
54 



cos 6(1, cos 6^ — cos 6g sin 6,f, sin 6^ 
cos 6(1, sin 6^ + cos 6g sin 6(f, cos 6^ 
sin 6(j, cos 6^ + cos 60 cos 6,/, sin 6^ 
sin 6(j, sin 6^ — cos 6g cos 6(f, cos 6^ . 



(31) 
(32) 
(33) 
(34) 



It can be verified using an algebraic calculation software that 9 and (/) given by equations 
(29) and (30) satisfy equations (20) and (21) with ci = 0, what means that the latter 
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equations describe correctly the precessional motion around an axis oriented arbitrarily in 
space. 

Before proceeding to the solution of equations (20)-(23) let us consider the elasticity 
of the Earth. This was first taken into account by Love [Tj and Larmor [2j because the 
predicted free precession period was of 305 days while the observed value was of around 
435 days. The correction due to Love introduces products of inertia C13 and C23 in Euler 
equations (8)-(9) for the Earth (see Kaula [16j) which become 

AlOi + (C - A)uJ2UJ3 - C23w| + CisWs = Ti (35) 

Alo2 - {C - A)ujiuj3 + Ci3u;| + C23UJ3 = T2 (36) 

with 

~ — 3G — 

- ^2 (38) 

where Re is the Earth mean radius and k2 is the so called Love number. Therefore we 
obtain result 

k2R%ujl\ . , k2Rluj'\\ 

k2R%uj'\\ / , k2R%uj'\\ 

By comparing this with equations (8)-(9) we see that the overall effect upon the last ones 
by the introduction of the elasticity of the Earth is to promote the transformation 

. , k2R%'Ui% / , , X 

A^A+ \ (41) 

The Love number k2 is left as one of the free parameters of this model. 

We solve numerically the set of four second-order differential equations given by (20)- 
(23), after performing the transformation given by equation (41), through the Runge- 
Kutta-Fehlberg (RKF) method. The time step upper limit is 1 solar day/1000 = 86.4 s; 
changing this limit to 1 solar day/10000 did not cause meaningful modifications in the 
results. We have 11 free parameters, i. e., the Love number k2, the momenta of inertia 
A* , C* of the internal ellipsoid, and the following that establish the initial conditions: the 
components c<j|*, u)2*, (^a* (recall that 0^3, is constant) of the angular velocity u* in the 
principal axes 1*2*3* of the internal ellipsoid, the Euler angles O12 = a, </'i2; V'12 that give 
the spatial orientation of the 1*2*3* frame with respect to the 123 frame, and the angular 
velocities (j) and (j)* . On the other hand, from equations (14) and (15) we have 



d = ±^vol + u:l-(tPsin^e (42) 

from which we can determine the initial value of 9 unless by its sign. In an analogous way 
we have 

r = ±Jujl^ + uj^i - </)*2 sin2 9*. (43) 
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^From equations (42)- (43) we can also find bounds for tlie initial guesses of (j) and (/>* 



<<P< , (44) 



sm y sm ( 



<<P*< • (45) 



smf^ smf 



Besides, from energy conservation of the Earth-ellipsoid system, we have the following 
relation to be satisfied by the initial guesses of and uj2* 

+ ^.2 > I j^^g(2a) - 1] - A (^2 + _ ^ (Ci3u;i + (46) 

where a is the initial guess for the angle between the axes 3 and 3*, and loi and L02 are 
the initial components of the Earth angular velocity u in the axes 1 and 2, respectively. 

Hopfner |1H [T2] has managed to filter the Chandler wobble (CW) and the annual 
wobble (AW) as major components from the data related to the Earth polar motion of 
its rotation axis around the polar principal axis. These data have being accumulated for 
more than 100 years and a compilation of them is published periodically by Gross [17j 
from NASA-JPL. We use a relation between the polar motion, described as the two angles 
PMX and PMY , and the orientation of the Earth rotation axis, described by the angles 
9^ and 4)^^ (see Fig. 5). As usual, the 1-axis is towards the Greenwich meridian and the 
2-axis is towards the 90°E meridian. The angles O^^j and t/)^ can be calculated from 



tan = ± V tan2 PMX + tan2 PMY (47) 
tanPMF , 

tan 6^ = (48) 

tan PMX 

Therefore by knowing PMX and PMY at the start of a Chandler wobble, taken from 
Hopfner [11] calculations, and by knowing the value of the length-of-day (LCD) at that 
moment, taken from Gross [IT] data, we can calculate the initial values of wi and u;2, which 
are necessary in equations (42), (44) and (46), and the value of ^3, which is considered as 
constant during the motion. Furthermore, we can calculate then the initial components 
of the Earth angular momentum L in the direction of the axes 123 as 

Li = Auji + Ci^uj^ (49) 
L2 = Auj2 + C23W3 (50) 
L3 = Ccja + Ci3tJi + C23W2- (51) 

Once we have guessed the initial values of uj\* , UJ2* , ^3, , we can also calculate the initial 
values of the components of L* in the axes 1*2*3* by 

LI, = A*Lol, (52) 
L*2* = A*L0*2^ (53) 
LI, = C*uj*,,. (54) 
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With these quantities available we can calculate the components of the total angular 
momentum in the axes 123 



Lti = Li + Ll = Li + L*, (cos (f)i2 cos ^"12 — cos a sin (j!)i2 sin ■012 ) 

— L2* (cos 012 sin ijji2 + cos a sin ^12 cos ■012) + -^3* sin a sin (^12 (55) 
Lt2 = L2 + L2 = L2 + L^. (sin cos -012 + cosq;cos0i2 sin ■012) 

+ L2* (— sin 012 sin ^^12 + cos a cos 0i2 cos -012) — L3* sin a cos 0i2 (56) 
Lts = L3 + L3 = L3 + Ll» sin a sin ■0i2 + L2* sin a cos -012 + L^* cos a. (57) 

Here, a = 612, 0i2 and ^^12 are the initial guesses for the Euler angles that describe the 
1*2*3* reference frame with respect to the 123 frame. Since we define the Z-axis of the 
fixed reference frame in the same direction of (see Fig. 6), the initial value of the Euler 
angle 6, that is necessary for the solution of equations (20)-(23), is obtained from 

n • U3 Lj3 

cos U = — — — = (58) 

* y^ti + ^t2 + ^3 

where U3 is a unit vector in the direction of the 3-axis. In order to calculate the initial 
value of we construct the Y orthogonal to the L*L-plane (Fig. 6), taking the 

initial values of such vectors, and the X-axis orthogonal to the YZ-plane, as usual. Then 

U3Y Lt2Li - LtlL2 , . 

cos = = — , . (59) 

U3XY l-L X Lt| sm 

where u^y and u^xv are the projections of U3 in the Y-axis and in the XY-plane, respec- 
tively. In order to determine completely we have also to express 

. , U3X Lt2{Lt2Ls — Lt3L2) — Lti{Lt3Li — LtiL^) , . 

sm0 = = ^ — -. (60) 

U3XY Lt\L xLt\sme 

In an analogous way we have the projections of in the axes 1*2*3* 

Lti* = Ll* + Li* = L^* -I- Li (cos 012 cos ^12 — cos a sin 012 sin -012) 

-I- L2(sin 012 cos ^12 + cos a cos 0i2 sin ^12) + L3 sin a sin ^12 (61) 
Lt2* = L2* + L2* = L2* — Li (cos 012 sin ^12 + cos a sin 012 cos -012) 

-h i>2(— sin 012 sin il)i2 + cos a cos 0i2 cos ^12) + sin a cos ^^12 (62) 
Lt3* = L3* -|- L3* = L3* -|- L\ sin a sin 0i2 — L2 sin a cos 0i2 + cos a. (63) 

and the initial values for the Euler angles 9* and 0* are obtained from 

• U3 Lt3* 

cost' = — j-^ = I = (64) 

* V ^a* + ^h* + ^3* 

,* Lti* L2* — Lt2* LI* 

cos = = — p— r i ■ a* (6^) 

u^xY |L*xLt|sm0* 
• A* — _ ^t2*{Lt3*L2* — Lt2*L%*) — Lti*{Lti*L\» — Lt3*L\*) 

^^'^^ ~ ulxY~ Lt|L* X Lt|sin6'* ' ^ ' 
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While integrating equations (20)- (23) we can calculate the CW motion through (see 
Fig. 5) 



tanPMX = — (67) 

^3 

U!'? 

tan PMY = — 

with uJi and uj2 given by equations (14) and (15). The necessary values of ip are obtained 
step- by-step through numerical integration of equation (16). 



Ill Results and discussion 

With the aim of getting 10 unknown parameters of our model (all but /C2) we find that 
it is enough to establish the following conditions: 1) reproduce the Chandler component 
of the Earth rotation axis position predicted by Hopfner [11] at the start of a wobble; 2) 
minimize the length of the prograde (counter-clockwise) trajectory followed during that 
wobble motion. The second condition is necessary in order to avoid intrusive nutations that 
might emerge accompanying the precessional motion. This means that we minimize, for a 
given initial guess of k2, the sum of: a) the squared distance from the model calculed Earth 
rotation axis position (PMX , PMY) at the start of the wobble and the one predicted by 
Hopfner; b) the wobble trajectory length in the plane {PMX, PMY) obtained by numerical 
integration. For this minimization process we used the POWELL routine from the Numerical 
Recipes package (Press et al. [ISj). However the generalized simulated annealing (GSA) 
code (Mundim and Tsallis [19] , Dall'Igna Junior et al. [2^ ) has proved to be very effective 
in a first step run in order to prevent getting stuck in local minima. The 11th parameter, 
is then found by reproducing the Chandler component of the Earth rotation axis 
position at the end of the wobble, with reference to Hopfner's prediction, so that ^2 is 
responsible for the adjustment of a particular wobble period. We show in Figs. 7-9 some 
of our results for the CW in three different epochs, in comparison with the results due to 
Hopfner [11], which are represented by dots. The first one, for which Wg* is negative, has 
the typical quality of providing reasonable reproduction of the motion. The next Figures 
show the two exceptions to this rule, which tend to present almost circular trajectories 
instead of, in those epochs, the slightly elliptical ones due to Hopfner. This may be caused 
by the above mentioned criteria - mainly condition 2) - adopted in the fit process. We also 
show in Table I the corresponding values for the 11 parameters (with respective numerical 
uncertainties shown in Table II), while in Table HI we have the initial values of 9 and 
9* calculated from equations (42) and (43) with appropriate signs and the corresponding 
values of the torque strength ci. As expected - see comments after equation (7) - a 
turns out to be close to 90°. On the other hand, parameters Wg* (the main angular 
velocity component of the ellipsoid). A* , C* and ^2 present oscillations when examining 
the respective columns in the Table I, being more noticeable for the 44330-44755 interval. 
For /c2, since it is connected to the period and that particular wobble has a short one 
- see Table IV for a comparison between the presently calculated periods and the ones 
due to Hopfner which they shall reproduce - this is not surprising although undesirable. 
Hopfner [IIIIT2], like other authors (see e. g. Liu et al. [21], Wang [22]), found variable 
CW periods but this is not a clearly solved question since some investigators have another 
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opinion (see e. g. Vicente and Wilson |23j . Jochmann |24) . Liao and Zhou [25j . Guo et 
al. [26\). For the sake of comparison, values of k2 obtained by some other authors are 
0.284 (Kaula pj6]), 0.30088 (model 1066A of Gilbert and Dziewonski ^) which are not 
far from what we got with the exception of the 44330-44755 case. 

By comparing the values of C* with momenta of inertia of natural satellites present 
in the solar system we find that they range from about Miranda's (Uranus) to Charon's 
(Pluto). Since we are able only to calculate the DM ellipsoid principal momenta of inertia, 
which are related to the ellipsoid mass M* and their equatorial and polar radii a* and c* 
through A* = M*(a*^ + c*^)/5 and C* = 2M*a*^/5, we cannot establish the ellipsoid mass 
uniquely, and in consequence, the fraction of DM present in the Earth is not predicted. 
Notice that, for an oblate ellipsoid, we have c* < a* and, therefore, the relation C*/2 < 
A* < C* holds and it was required to be satisfied in our numerical procedure. Just to 
have an idea of the order of magnitude of M*, let us suppose that a* ^ 2350 km, i. e., 
the ellipsoid has the average radius of the outer core of the Earth. From the values of C* 
in Table I, this implies that the fraction of DM present in the Earth would lie in the range 
6 X 10^^ to 9 X 10~^. We emphasize that, according to the model, the Earth motion and 
its gravitational interaction with other bodies have already incorporated the presence of 
DM through what has been known effectively as "Earth mass" , so that the outcome of its 
presence might be more easily observed through the Earth wobble. 

Following the hypothesis of changing CW periods we present in Fig. 10 the corre- 
lation between period and average amplitude. We show Hopfner's results [TT] and ours. 
The last ones are fitted by the full line, which has the form: period (in days) = 440- 
3101. 47xexp[-0. 0347946 x amplitude (in mas)]. As a reference we also present the fits 
due to Liu et al [21] related to the epochs 1912-1928 and 1936-1948. The correlation 
clearly varies with time. 

One might argue that the parameter fluctuations in Table I indicate that this model is 
not sound. However those fluctuations may be caused by the model simplistic - although 
convenient at this stage - approach that we are considering the DM ellipsoid the only source 
responsible for the complexity of Chandler wobble, what obviously is not true. Therefore, 
the predicted A* and C* dark matter ellipsoid parameters should be considered only as 
upper limit values. 

Another point that could be raised against the model is that it is just emulating one or 
more of the attributed causes for CW - as listed in the Introduction - and, in consequence, 
the DM ellipsoid existence is not real. However, it can be shown that externally to the 
DM ellipsoid, whose rotation axis is almost orthogonal to the Earth's, it generates a non- 
isotropic, time dependent, gravitational fleld with two contributions, with respect to a 
U1U2U3 reference frame attached to the Earth. The flrst term is radial: 



Here, 9 is the co-latitude, (p is the longitude and ujt is the Earth angular velocity. Of 
course, ujxt is deflned except for a phase that would establish the position of the ellipsoid 
rotation axis with respect to the Greenwich meridian (for instance, in the direction of ui), 
at a given time. The second contribution is 



3G{C* -A*) r 
2^ I 



3 sin^ d{cos (f) sin uJTt + sin cos uoTt)"^ — 1 • 



(69) 



a.Q* — 



?,G{C* -A*) 
2?^ 



sin 20* [{cosO* sm(j)* cosLOTt — smO* sm.LOTt)'^! 
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— (cos 9* sine/)* siuLdTt + sin 9* cos W7-t)u2 + cos 9* cos <p*n^] (70) 



where 

sm9* = \J cos^ 9 + sin^ 9 [cos (f) cos — sin(/)sin wj't)^ (71) 

COS0* = sin0(cos(/)sinL<JTt + sin^cosaJT^)) (72) 

sm0(cos(/)cosLiJTi — siiK^isinwri) /„ x 

smcp = = (73) 

y cos^ + sin^ 9 (cos (/) cos uort — sin (/) sin ooTt)"^ 

coscp = . (74) 

w cos^ 9 + sin^ 0(cos (p cos uxt — sin (/) sin uJxtY 

Note that, in the Earth's reference frame, the DM elhpsoid gives a complete turn in a 
day. Supposing the DM elhpsoid is smaller than the Earth core, this gravitational field 
certainly enforces viscous flow in the outer core, which is in a fluid state. As a first 
consequence, it can provide the necessary energy to the magneto-hydrodynamic motion 
in the outer core, which is closely related to the generation of the geomagnetic field. This 
could give a plausible explanation to one of the basic problems of geophysics (see Stacey, 
ref. [I3]), namely, what is the source of this energy. However the numerical solution 
of this motion through the Navier-Stokes and energy conservation equations suffers lack 
of knowledge of the flowing material density dependence with respect to pressure and 
temperature. Secondly, it can generate heat in the core, that would contribute to decrease 
the shortfall of 0.7 TW to the energy necessary to maintain the adiabatic temperature 
gradient in the core, leaving more power to drive the geomagnetic field dynamo [13]. On 
the other hand. Mack et al. |28j have concluded that DM is unlikely to contribute not only 
to Earth's internal heat fiow but also to hot-Jupiter exoplanets. However they studied the 
contribution of DM self-annihilation to heat generation, what is a different idea from ours. 
Conversely, Adler [9] has also studied the contribution of self-annihilating and non-self- 
annihilating DM accretion to the internal heat of the Earth, Jovian planets and hot- Jupiter 
exoplanets. His conclusion is that this process is plausible provided efficient DM capture 
is occurring. In our description, it is feasible to imagine DM ellipsoids also present is such 
bodies, generating heat by enforcing internal viscous flow, like in Earth. Therefore this 
model can mean much more than mere emulation of known suggested causes for CW but 
it also touches open questions like geomagnetic field dynamo and heat generation in the 
Earth's outer core and in other planets. 

Rather than giving the final answer for the problem that has raised the attention of geo- 
physicists and astronomers for several decades, we expect that this calculation may open 
room for a new and potentially important key component in CW which was completely 
unsuspected untill now, to be considered in more sophisticated models. Future devel- 
opment will decide about the real significance of this approach that might have relevant 
consequences in our comprehension about dark matter. At least, it has been demonstrated 
that two interacting ellipsoids can have Chandler-like wobble. 
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FIGURE CAPTIONS 



Figure 1 - Two ellipsoids tilted by an angle a around the common equatorial principal 
axes 1 and 1*. The spherical coordinates (r, 9, (p) that give the position of an infinitesimal 
mass dm with respect to the 1*2*3* reference frame are also shown. 

Figure 2 - Definition of the Euler angles 9, (p, ip. 

Figure 3 - Space and body cones generated by the Eulerian (free) precession. The axes 
1, 2 (equatorial principal axes) and 3 (polar principal axis) are fixed in the oblate ellipsoid 
and follow its motion. The angular momentum L defines a fixed direction in space while 
the angular velocity u gives the instantaneous rotation axis. 

Figure 4 - Free precession around an arbitrary axis z. The xyz-frame is oriented with 
respect to the XYZ-frame, which is fixed in space, through the Euler angles 5e, Stf,, 5^. 

Figure 5 - Definition of the angles PMX and PMY which describe the Chandler wobble 
motion. The origin of the 123-frame is at the Earth center. Axes 1 and 2 are towards 
the Greenwich meridian and 90°E meridian, respectively. The 3-axis is the polar principal 
axis. The angular velocity uj is also shown as well as its spherical coordinates 9^ and 0^ . 

Figure 6 - Definition of the frame XYZ fixed in space which is convenient to describe 
both interacting ellipsoids motion. The Z-axis coincides with the direction of the total 
angular momentum L^. The XZ-plane is the same plane formed by the L (Earth) and L* 
(ellipsoid) angular momenta. The unit vectors U3 and U3, which are in the direction of 
the polar principal axes 3 and 3*, are also shown as well as their spherical coordinates 0, 
(t>s, 9*, </.*. 

Figure 7 - Calculated (line) and predicted (dots) Chandler wobble prograde (counter- 
clockwise) motion in the interval 44330-44755 (in modified Julian date - MJD) or 01/APR/ 
1980-31/MAY/1981. The asterisk shows the position of the motion start point. PMX and 
PMY are given in milli-arcsec (mas). 

Figure 8 - Same as in Fig. 7, for the interval 44970-45403 (MJD) or Ol/JAN/1982- 
lO/MAR/1983. 

Figure 9 - Same as in Fig. 7, for the interval 50170-50604 (MJD) or 28/MAR/1996- 
05/JUN/1997. 

Figure 10 - Relation between CW period (in days) and average amplitude (in mas). 
Circles represent results due to Hopfner [TT| while triangles show our results. The full line 
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fits the last points. As a reference, daslied and dotted lines represent fits for the epochs 
1912-1928 and 1936-1948 respectively. 

TABLE CAPTIONS 

Table I - Results for the 11 unknown parameters in several epochs: initial values for 
the angular velocities </> and (p* (in rad/s), for the components uJ2*, (in rad/s) 
of angular velocity uj*, for the Euler angles a = 612, (j)i2, ipi2 (in degrees); values of the 
internal ellipsoid principal momenta of inertia A*, C* (in kg • m^) and Love number k2. 

Table II - Numerical uncertainty estimates (in %) of the parameters shown in Table I. 

Table III - Calculated initial values of 9 and 9* (in rad/s), with respective appropriate 
signs, and torque strength ci (in N • m^). 

Table IV - CW periods resulting from the model in comparison with the values from 
Hopfner [11] (in days). 
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Table III 



interval 


a 




e 




Cl 


44330-44755 


2.426E- 


-11 


4.141E- 


-10 


9.226E-h23 


44970-45403 


-2.684E- 


-11 


1.419E- 


-07 


1.367E-F22 


45700-46135 


-4.176E- 


-11 


7.080E- 


-07 


6.790E-h20 


46340-46773 


3.218E- 


-11 


3.487E- 


-06 


5.006E-h20 


46977-47410 


-2.231E- 


-11 


6.839E- 


-07 


6.349E-F20 


47640-48075 


3.367E- 


-11 


2.891E-06 


3.922E-F20 


48260-48696 


-4.771E- 


-12 


-4.521E- 


-09 


2.607E-F21 


48900-49336 


-9.816E- 


-12 


1.147E- 


-08 


2.517E+21 


49630-50066 


-5.563E- 


-11 


3.185E- 


-10 


3.331E-h20 


50170-50604 


-4.110E- 


-11 


4.946E- 


-10 


1.900E-F21 



Table IV 



interval calculated period Hopfner's period 



44330-44755 


424.578 


424.76 


44970-45403 


432.454 


432.45 


45700-46135 


434.677 


434.68 


46340-46773 


433.366 


433.37 


46977-47410 


433.331 


433.33 


47640-48075 


435.123 


435.12 


48260-48696 


436.414 


436.41 


48900-49336 


436.943 


436.93 


49630-50066 


436.924 


436.92 


50170-50604 


433.920 


433.92 
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Table I 



interval 









LO2* 






(f>12 


■012 


a 


A* 


c* 


k2 


44330-44755 


1.386E-7 


3.087E-20 


9.907E-12 


4.140E-10 


-1.592E-4 


358.849 


353.219 


90.001 


6.334E-F32 


1.264E-F33 


0.2672 


44970-45403 


-2.704E-5 


-9.685E-7 


-9.789E-7 


-2.088E- 


-9 


4.952E-4 


173.698 


337.477 


89.911 


2.399E-h31 


3.334E-F31 


0.2801 


45700-46135 


5.711E-5 


-4.203E-6 


-4.262E-6 


-2.128E- 


-9 


6.853E-4 


359.988 


359.999 


89.847 


1.079E-F31 


1.125E-h31 


0.2835 


46340-46773 


-7.419E-5 


1.823E-6 


3.935E-6 


-2.410E- 


-9 


1.089E^3 


187.700 


350.276 


89.849 


7.429E-F30 


7.772E-F30 


0.2815 


46977-47410 


7.205E-5 


-4.490E-6 


-4.542E-6 


-2.128E- 


-9 


7.455E-4 


359.988 


359.999 


89.857 


1.168E-I-31 


1.211E-F31 


0.2815 


47640-48075 


-3.769E-5 


2.098E-6 


3.572E-6 


-2.410E- 


-9 


7.727E-4 


109.191 


353.481 


89.832 


7.280E-h30 


7.548E-h30 


0.2842 


48260-48696 


5.317E-6 


-1.861E-7 


-1.862E-7 


2.128E-- 


-9 


3.338E-3 


222.469 


359.999 


90.047 


1.947E+31 


2.125E+31 


0.2861 


48900-49336 


-2.505E-5 


-1.142E-6 


-1.142E-6 


-2.088E- 


-9 


5.632&-4 


346.867 


332.210 


89.841 


1.731E-F31 


1.903E-F31 


0.2869 


49630-50066 


4.133E-5 


-6.759E-6 


-6.759E-6 


-2.128E- 


-9 


7.390E-4 


359.988 


359.999 


89.873 


9.482E-I-30 


9.710E-F30 


0.2869 


50170-50604 


-3.317E-5 


-6.301E-7 


-6.301E-7 


-2.088E- 


-9 


5.230E-4 


247.397 


333.076 


89.877 


8.287E-h30 


9.587E-h30 


0.2824 



Table II 



interval 







A4>* 
















A</.i2 




Al/'12 




Aq 




AA* 




AC* 




Ak2 


4> 




















't>l2 




^"12 




a 




A* 




C* 




k2 



44330-44755 


0.10 


4.6E-6 


1.3E-7 


2.3E-9 


71 


1.6E^3 


1.6E-2 


6.6E-4 


l.lE-3 


129 


l.lE-4 


44970-45403 


1.8E-6 


0.95 


0.77 


3.0E-5 


2.5E-7 


2.2 


0.37 


5.0E-2 


4.1E-6 


28 


l.lE-3 


45700-46135 


9.6E-6 


1.5 


0.16 


1.8E-6 


1.6E-6 


4.5E-8 


1.9E-10 


9.1E-6 


4.3 


1.8 


1.2E-5 


46340-46773 


3.0E-4 


0.44 


0.11 


3.7E-5 


13 


16 


2.8E-2 


3.1E-4 


4.6 


7.0 


6.6E-5 


46977-47410 


1.2E-6 


1.3 


0.64 


2.7E-5 


6.2E-7 


7.2E-8 


<1E-13 


7.9E-6 


3.8 


3.5 


5.8E-5 


47640-48075 


3.6E-4 


3.8 


0.26 


2.0E-5 


1.6E-5 


6.2 


3.5E-2 


7.4E-5 


4.1 


1.3 


1.2E-4 


48260-48696 


3.3E-5 


2.8E-2 


1.7E-3 


3.4E-5 


38 


26 


7.7E-10 


7.0E-6 


14 


22 


4.9E-5 


48900-49336 


8.3E 5 


1.1 


1.9 


5.0E-5 


4.6 


0.34 


1.4 


1.5E-5 


9.9 


13 


1.6E-4 


49630-50066 


<1E-13 


2.6 


6.8E-2 


2.6E-7 


1.6E-3 


7.9E-14 


<1E-13 


3.2E-9 


2.7 


0.38 


2.5E-5 


50170-50604 


2.7E-5 


2.3 


l.OE-6 


3.0E-6 


<1E-13 


l.lE-2 


2.2E-2 


2.7E-9 


18 


11 


1.3E-4 
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